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Abstract In this paper we deal with Mellin convolution of generalized Gamma densities which 
leads to integrals of modified Bessel functions of the second kind. Such convolutions allow us to ex- 
plicitly write the solutions of the time-fractional diffusion equations involving the adjoint operators of 
a square Bessel process and a Bessel process. 
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1 Introduction and main result 

In the last years, the analysis of the compositions of processes and the corresponding governing equa- 
tions has received the attention of many researchers. Many of them arc interested in compositions 
involving subordinators, in other words, subordinated processes Y(T(t)), t > (according to [9j) where 
T(t), t > is a random time with non- negative, independent and homogeneous increments (see Q). 
If the random time is a (symmetric or totally skewed) stable process we have results which are strictly 
related to the Bochner's subordination and the p.d.e.'s connections have been investigated, e.g., in 
@; S B EEl • If the random time is an inverse stable subordinator we shall refer to the govern- 
ing equation of Y(T(t)) as a fractional equation considering that a fractional time-derivative must be 
taken into account. In the literature, several authors have studied the solutions to space-time fractional 
equations. In the papers by Wyss (30| . Schneider and Wyss [2{J, the authors present solutions of the 
fractional diffusion equation d^T = d^T in terms of Fox's functions (see Section [2]). In the works by 
Mainardi et al., see e.g. [l?]; EH the authors have shown that the solutions to space-time fractional 
equation x DqU = tD*u can be represented by means of Mellin-Barnes integral representations (or 
Fox's functions) and M- Wright functions (see e.g. Kilbas et al. [13|)- The fractional Cauchy problem 
Dfu — Lu has been thoroughly studied by yet other authors and several representations of the solu- 
tions have been carried out, but an explicit form of the solutions has never been obtained. Nigmatullin 



25 1 gave a physical derivation when L is the generator of some continuous Markov process. Zaslavsk 



sky 

m 



3l| introduced the space-time fractional kinetic equation for Hamiltonian chaos. Kochubei 14, 
first introduced a mathematical approach while Baeumer and Meerschaert [l[ established the connec- 
tions between fractional problem and subordination by means of inverse stable subordinator when L 
is an infinitely divisible generator on a finite dimensional vector space. In particular, if dtp = Lp is the 
governing equation of X(i), then under certain conditions, q = Lq+5(x)t~^ /T(l — j5) is the equation 
governing the process X(V t ) where V* is the inverse or hitting time process to the /3-stable subordinator, 
(3 <= (0, 1). Orsingher and Beghin [2(1 [27| found explicit representations of the solutions to d"u = \ 2 d^u 
only in some particlular cases: v = (1/2)", n e N and v — 1/3,2/3,4/3. Also they represented the 
solutions to the fractional telegraph equations in terms of stable densities, see [3j;[26|. In general, the 
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solutions to fractional equations represent the probability densities of certain subordinated processes 
obtained by using a time clock (in the following we will refer to it as L v t ) which is an inverse stable sub- 
ordinator (see Section [4]). For a short review on this field, see also Nane [13] and the references therein. 

We will present the role of the Mellin convolution formula in finding solutions of fractional diffusion 
equations. In particular, our result allows us to write the distribution of both stable subordinator 
and its inverse process whose governing equations are respectively space-fractional or time-fractional 
equations. This result turns out to be useful for representing the solutions to the following fractional 
diffusion equation 

BlnY = Q^uY (1.1) 

where u 7 /' 1 — ul'^(x,t), x > 0, t > 0, is the Riemann-Liouville fractional derivative, v G (0,1] 
and G 1:l i is an operator to be defined below (see formula (|3.4jl ). We present, for v = l/(2n + 1), 
n G MU {0}, the explicit solutions to in terms of integrals of modified Bessel functions of the 

second kind (K u ) whereas, for v G (0, 1], we obtain the solutions to in terms of Fox's functions. 
After some preliminaries in Section[2l in Section[3]we recall the generalized Gamma density starting 
from which we define the distribution of the (generalized Gamma) process Gj' 11 and the distribution 
e 7 of the process E^ 4i . The latter can be seen as the reciprocal Gamma process, indeed E^'^ = 
or in a more striking interpretation, as the hitting time process for which (i?7' M < x) = (G 7,M > t). We 
shall refer to i?7' M as the reciprocal or equivalently the inverse process of G]'^. It must be noticed that 
e Ji — 9/7 7 because G^ 1 '* 1 = 1/Gj'^. Furthermore, we introduce the most important tool we deal with 
in this paper, the Mellin convolutions g^* n (see formula (|3.14[) ) and e-ff 1 (see formula p,13|) ) where 
e*^ stands for e^'*". In Section 2] we draw some useful transforms of the distribution h v of the stable 
subordinator t% and the distribution l v of the inverse process L" . Similar calculations can be found 
in the paper by Schneider and Wyss [29j. The inverse (or hitting time) process is defined once again 
from the fact that (L v t < x) — (f£ > t) (see also [J 111)- In Section [5] we present our main contribution. 
We show that the following representations hold true: 

h v {x,t) = ep n (x,(p n+1 (t)), x > 0, t > 0, v = 1/(71 + 1), neN 

and 

l u (x,t) = g£ +1) '* n (x,i> n+1 (t)), x>0,t>0, i/ = l/(n+l), neR 

where /2 = (fi\, . . . , fj, n ), fj,j = j v , j = 1, 2, . . . , n, v = l/(n+l), n G N and the time-stetching functions 
are given by <p m (s) — (s/m) m and ip m (s) — rns 1 /™, s G (0,oo), m G N, ip = f^ 1 - 
The discussion made so far allows us to introduce the result stated in Theorem [TJ For v = l/(n + 1), 
«eNU {0}, the solutions to (jl.ip can be written as follows 

/>oo 

fiZ'^M) - / 9lt*,s lh ) gy" MX/u ~ X \s^ 1/v {t)) ds, x G (0,oo), * > 
Jo 

where, for n G 2N U {0}, we have 

1 / x \ 1 ~ f°° f°° 

g X/u,*{X/u-X)^^ = __ *» ^ ^ Q i _ J ,( il . (fll )...Q i _ :: ( Sn _ 1>t )d fll ...d fln _ 1 

and 

Qi^l {x, t) = ifw U^x/tf/A , x > 0, t > 0. 

As a direct consequence of this result we obtain iij'^ = g 7 , for v = 1, where gl{x, t) — g^(x, i 1 / 7 ) and 
the governing equation writes 
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Furthermore, for 7 = 1,2 and v e (0, 1] we obtain 



D »*=(*|l-<"- 2 >J-)*' 



x > 0, t > 0, fi > (1.2) 



and 

D «" = MS-S^)^ *>M>o,„>o. (1.3) 

Equation p. 31) represents a fractional diffusion around spherical objects and thus, the solutions we 
deal with obey radial diffusion equations. 



2 Preliminaries 



The H functions were introduced by Fox [l0( in 1996 as a very general class of functions. For our 
purpose, the Fox's H functions will be introduced as the class of functions uniquely identified by their 
Mellin transforms. A function / for which the following Mellin transform exists 



M[f(-)](r)) 



x"f(x 



dx 



can be written in terms of H functions by observing that 

dx 



n Trm.n 



5R{r?} > 



where 



UT=i T (bj + vfr) Uti r(i - oi - mi) 

nLm+i r ( i - h - vPj) nr=n+i r K + ^ 



(2.1) 



(2.2) 



The inverse Mellin transform is defined as 



2ni 



6+ioo 



f(x) = 7T- / M\f(>)](v)x-"dn 



9—ioo 



at all points x where / is continuous and for some real 0. Thus, according to a standard notation, the 
Fox H function is defined as follows 



Tjm,n 



1 

2iri 



MZfWx-^dr, 



P{V) 



where P(X>) is a suitable path in the complex plane C depending on the fundamental strip (2?) such 
that the integral (|2.1[) converges. For an extensive discussion on this function see Fox Mathai and 
Saxena [20(. The Mellin convolution formula 



ds 

fi*fi(x)= I fi(x/s)f 2 (s)—, x>0 



turns out to be very useful later on. Formula (|2.3p is a convolution in the sense that 

M [fi*h{-)\ (v) =M[fi(-)} ( V ) x M [h(-)\ fa). 
Throughout the paper we will consider the integral 



fi°f2(x,t)= f 1 (x,s)f 2 (s,t)ds 
Jo 



(2.3) 



(2.4) 



(2.5) 
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(for some well-defined fx, $2) which is not, in general, a Mellin convolution. We recall the following 
connections between Mellin transform and both integer and fractional order derivatives. In particular, 
we consider a rapidly decreasing function / : [0, 00) i— > [0, 00), if there exists a € K such that 



lim x a - K - x -^-rf(x) = 0, fc = 0,1,..., n-l, n G N, x € R+ 



then we have 



M 



JTL 
— /0 

dx nJy 



(*?) =(-l) n p^l^ I/O)] (»?-n) 
r(ry - 71) 



and, for < a < 1 



r(7?- a) 



(2.6) 



(2.7) 



(see Kilbas et al. [X3j | : Samko et al. [28[ for details). The fractional derivative appearing in (|2.7j) must 
be understood as follows 



"7-/0*0 = TT7 r / (x - s)"~ Q - 1 ^(s) da, n - 1< a < n 

dx a w r (71 - a) 7 v ; ds" w 



(2.8) 



that is the Dzerbayshan-Caputo sense. We also deal with the Riemann-Liouville fractional derivative 



D a f = 

X J 



r (n — a) dx n Jq 



(x - s y 



1 f(s) ds, n — 1 < a < n 



and the fact that 



k=0 



x=0+ 



r(fe-a + l)" 



n — 1 < a < n, 



(2.9) 



(2.10) 



see Gorenflo and Mainardi [ll| and Kilbas et al. [13|. We refer to Kilbas et al. [13|; Samko et al. [26 
for a close examination of the fractional derivatives (12.81) and 



3 Mellin convolution of generalized Gamma densities 

In this section we introduce and study the Mellin convolution of generalized gamma densities. In the 
literature, it is well-known that generalized Gamma r.v. possess density law given by 

Ql{z) = jjr^exp{-z^}, z>0, 7 > 0, /i > 0. 



Our discussion here concerns the function 
Let us introduce the convolution 



exp 



x > 0, t > 0, 7 ^ 0, n > 0. 



- ^ 2 (x, t) = / (z, a^J* (s, t)d* = sign( 7l72 ) 



1 Jo 



ds 



for which we have (see formula (|2.4p ) 

■m [<$ *<$(■,*)] fa) = m ki(->* 1/2 )l to) x ^ to->* 1/2 )l (»?) 



(3.1) 



(3.2) 



(3.3) 
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as a straightforward calculation shows. We now introduce the generalized Gamma process (GGP in 
short). Roughly speaking, the function (|3.1[) can be viewed as the distribution of a GGP {G]'^, t > 0} 



in the sense that Vi the distribution of the r.v. G] is the generalized Gamma distribution (|3.ip . 
Thus, we make some abuse of language by considering a process without its covariance structure. 
In the literature there are several non-equivalent definitions of the distribution on K™ of Gamma 
distributions, see e.g. Kotz et al. [3] for a comprehensive discussion. In Section [5] (Corollary [l| we 
will show that the distribution (|3.1I) satisfies the p.d.e. 



d_ 7 

at 9 » 



where 



2-7 



dx 



dt 

d_ 

dx 



(7M - 1) 



x > 0, t > 



„l-7 



dx" 



/, x > 0, t > 



(3.4) 



and 7 ^ 0, / £ D(Qy tfl ). For 7=1, equation (|3 . 1 1) becomes the distribution of a 2/j-dimensional 
squared Bessel process {BESSQ^$ , t > 0} and, for 7 = 2 we obtain the distribution of a 2/i- 
dimensional Bessel process {BESftf, t > 0}, both starting fr om zero. Some interesting distributions 



can be realized through Mellin convolution of distribution g^. Indeed, after some algebra we arrive at 

7 2,7/^1-1^7^2 



and 



<(x,t) 



h*9^ (x,t) 



ry 3,7^2-1^7^1 



x > 0, t > 0, 7 > 



x > 0, t > 0, 7 > 



(3.5) 



(3.6) 



B(nx,fj, 2 ) (n + x7)Mi+M2^ 

where B(-, •) is the Beta function (see e.g. Gradshteyn and Ryzhik 12J, formula 8.384]). Moreover, in 
light of the Mellin convolution formula (|2.4[) . the following holds true 

M [gl *g£(;tj\ („) - M [9^*9^ (-,«)] (ry). 

A further distribution arising from convolution can be presented. In particular, for 7^0, we have 



9l^9l 2 {x,t) 



2| 7 | (^/r)— 
xr(/Ji)r(/j 2 ) 



- K 




2; > 0, i > 



(3.7) 



which proves to be very useful further on. The function K v appearing in (|3.7[) is the modified Bessel 
function of imaginary argument (see e.g 12J, formula 8.432]). For the sake of completeness we have 
writen the following Mellin transforms: 



3=1 

7 



and 



M[gZ(x,.)] ( V ) = 
Formula (|3.8p suggests that 



r(/i) 
rf M - 



fTJ-l 



r^-l 



t > 0, 5f?{r;} > 1 - 7M , 7 ^ 0, 



> 0, 5iM > 7/1, 7 ^ 0. 



(3.8) 



X [9l\ ■)] fa) = X [<?£(z 1/2 , •)] fa) x M [gllix 1 ' 2 , •)] (77). 

For the one-dimensional GGP we are able to define the inverse generalized Gamma process {E^ 
t > 0} (IGGP in short) by means of the following relation 

PriE]^ <x} = Pr{Gl<» > t}. 
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The density law = eZ(x, t) of the IGGP can be carried out by observing that 

d 



eUx.t) = Pr{E^ e dx}/dx 



and, making use of the Mellin transform, we obtain 



dx 



£f^(s, x) ds, x > 0, t > 



(3.9) 



M[el(-,t)] (r,) = 



M 



(r?) ds, ^i{r]} < 1 

/OO 
[#(*,•)] (»7-l)da 



= [by 



2-1 



rfa) 



The derivative under the integral sign in (|3.9|) is allowed from the fact that Si (s) = -^g^(s,x) € L 1 ! 
as a function of s. From (12.21) and the fact that 



H p,q 



H m,n 

p,g 



for all c > (see Mathai and Saxena [20j]), we have that 



n 


(M) " 




OM) _ 



(3.10) 



(3.11) 



By observing that M [e^(-,t)] (1) = 1, we immediately verify that (|3.11j) integrates to unity. The 
density law g 1 ^ can be expressed in terms of H functions as well, therefore we have 



x' 1 


0*,o) " 


~V 


(M) 



, x > 0, t > 0, 7 > 0. 



(3.12) 



In view of p. lip and (|3.12j) we can argue that 



Law s~i— 



1/G7 ,M , t > 0, 7 > 0, pL > o 



and e^(x,t) = gp(x,t), j > 0, x > 0, t > 0. 

Remark 1. We notice that the inverse process {E^' 1 ^ 2 , i > 0} can be written as 

£ t 14/2 = inf{s; B(s) = V2t} 
where B is a standard Brownian motion. Thus, E 1 ' 1 / 2 can be interpreted as the first-passage time of 



a standard Brownian motion through the level y 2t. 

In what follows we will consider the Mellin convolution e*[ l (x, t) = e Ml * . . . *e Mll (x, t) (see formulae 
(|2.4I) and p.2p ) where /2 = (/ii, . . . , /z„), /Xj > 0, j = 1,2, ...,n and, for the sake of simplicity, 
efj,(x,t) — e^(x,t). For the density law e*p{x,t), x > 0, t > we have 



r ( Mj + i-v) 



x [er (•.*)] (»?) = K-(-:* i/n )j w = f- 1 n x 

with 3?{?7} < 1. Furthermore, for the Mellin convolution g^'* n (x,t) = ffX*, • ■ • > *<7u„ O 2 ^ we have 



(3.13) 



n i 11 T (—+t*j) 

^br"(- ; *)]fa)=n^ ko.* i/n ) fa)=*"- i n 7 



i=i 



i=i 



rfa,) 



(3.14) 



with > 1 — miiijl/ij} 
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Lemma 1. The functions g 1 ^ and are commutative under convolution. 

Proof. Consider the Mellin convolution (|3.13l) . Let e». be the distribution of the process X a \ then 
formula (|3.13|) means that 



7,-1 



E {X 171 pT 2 (. . . X a " (t) . . =E^X ai (t^ n )X a2 (t 1 ^ 1 ) ■ • ■ X a " 

for all possible permutations of {o~j}, j = 1,2, ... ,n. The same result can be shown for e^. . Suppose 
now that the process X aj possesses distribution g 1 . , from (|3. 14|) we obtain the claimed result. □ 



4 Stable subordinators 

The ^-stable subordinators {f t , t > 0}, v G (0, 1), are defined as non-decreasing, (totally) positively 
skewed, Levy processes with Laplace transform 



and characteristic function 



where 



Sexp{-Af t M } = exp{-a l/ }, <> 0, A>0 
Sexp{z^ (l/) }=exp{-^(0}, 

du 



(4.1) 
(4.2) 



MO 



V 



„ 1 ' ; r(i-i/)w+i 



(see Bertoin Zolotarev (32j). After some algebra we get 

MO =<r\tf (l - * sgn(£) tan (^)) = , 



7T^ £ 



For the density law of the ^-stable subordinator {f t , t > 0}, say h v — h v (x,t), x > 0, t > we have 
the f-Mellin transforms 



M 



and 



2 |£| 

MA,-) (, 7 ) = A-^r(7 7 ) 



MM (v) = \Z\- n "e W ii^-^}r(ri) 



M 



(4.3) 



(4.4) 



where h u (^, t) = T [h v (-, t)] (£) is the Fourier transform appearing in (|4.2[) and h v (X, t) — C [h w (-, t)] (A) 
is the Laplace transform (|4.1j) . By inverting (|4.3|) we obtain the Mellin transform with respect to t of 
the density h v which reads 



_r(r?)T(l - rjv) f e iZL ^ 
T(rj)r(l- W ) 



(4.5) 



27T X 

r(r?)r(i-^) 



i — r\v) ( ( 7r -i r 7T ~n 
jcxp (-» o + w™?*' j + exp jz - - zttt^ j j 

, 2 > 0, i/e (0,1) 



■sin , T //,/ = I^-x""- 1 
rra: 1 -^ T(rp) 



where 5i{?7^} € (0,1). Formula (|4.5I) can be also obtained by inverting (|4.4[) . We are also able to 
evaluate the Mellin transform with respect to x of the density law h v . From (|4.3[) and the fact that 



J°° x v-i e -i^ dx= M. t where ( ± ^ = |^ ex p|±i^|||, us (0,1 



(4.6) 
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we obtain 



MM-.*)] (*?) = 



2 
2tt 



>7 ,,-**»(-£) 



2tw 



} 



1-7? 



v ) vY (1 — if) 



, km g (o, l), t > o. 



(4.7) 



The inversion of Fourier and Laplace transforms by making use of Mellin transform has been also 
treated by Schneider and Wyss [29j. 

We investigate the relationship between stable subordinators and their inverse processes. For a 
^-stable subordinator {f t , t > 0} and an inverse process {Lt"\ t > 0} (ISP in short) such that 

Pr{L^ ] < x} = Pr{f M > t} 
we have the following relationship between density laws 



l u (x,t) = Pr{L [ t v) G dzj/diz; 



<9x 



h v {s,x)ds, x > 0, t > 0. 



(4.8) 



We observe that -J^h u (s,x) exists and there exists £(s) € L 1 (]R + ) such that S 2 (s) = -^h u {s,x) = 
const ■ D^h u (s,x) < C( s )- The function /i„ is the distribution of a totally skewed stable process, 
thus /i„(x), x G R™ belongs to the space of functions in -D((— A) y / 2 ), see Samko et al. [28J. Thus, the 
integral in (|4.8j) converges. The density law (|4.8|) can be written in terms of Fox functions by observing 
that 



M [l v (;t)} fa) = / M 



(?7) 



/oc 
M[^(s,-)] (77-1; 



ds 



rfa) 



r (77^ — u + 1 

Thus, by direct inspection of (|2.2p . we recognize that 



i r (r)v - v) 
t"^-^, SR{ry> < l/v, t > 0. 



(l-I/,!/) 

(0,1) 



x > 0, t > 0, z/(0, 1). 



Density (|4.10[) integrates to unity, indeed M. [l v (-,t)] (1) = 1. The i-Laplace transform 

C[l v (x,-)](X) = \"- 1 exp{-x\ 1 '}, A > 0, v G (0,1) 
comes directly from the fact that 

J™ e- xt M [!„(•, t)} ( V ) dt = = J™ x^C [l v (x, •)] (A) dx. 



(4.9) 



(4.10) 



(4.11) 
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From (|4.11|) we retrieve the well-known fact that C [l u (-,t)] (A) = E v (—Xt u ) (see also Bondesson et al. 
[5]) where Ep is the Mittag-Leffler function which can be also written as 

EJ-\t v ) = - [°° exp {-^'"txX f lsm7TV dx t > 0, A > 0. (4.12) 

V TT J I J 1 + 2x v COS TtV + X 2u V ' 

The distribution l v satisfies the fractional equation -§^l v = —-§^vi x > 0, t > subject to l u (x, 0) = 
6(x) where the fractional derivative must be understood in the Dzerbayshan-Caputo sense (formula 
(|2.8[) ). The governing equation of l v can be also presented by considering the Riemann-Liouville 
derivative (|2.9I) and the relation (12.10)) (see e.g. Baeumer and Meerschaert Meerschaert and Scheffler 
[2~0 ] ; Baeumer et al. @ ) . It is well-known that the ratio involving two independent stable subordinator 
{ if t , t > 0} and { 2Tt V \ t > 0} has a distribution, Vi, given by 

r( W ) = Pr{ x t\ v) I 2 t\ v) G dw}/dw = - — ^\ 2v , w > 0, t > 0. (4.13) 

TT 1 + 2w v COS 7T!/ + W ZV 

Here we study the ratio of two independent inverse stable processes t > 0} and { 2 Ll u) , t > 0} 

by evaluating its Mellin transform as follows 

e{ x L^/ 2 L^Y~ 1 =M U;t)] (r,) x M U;t)] (2-rj) = ^"^"^ (4-14) 

with dt{r]} G (0, 1). By inverting (|4.14[) we obtain 



1 sin vtt 1 



. , sin VTT — 11VTT „ , , 

fcfa;) = - = — / : — x^drj (4.15) 

vtt 1 + 2x cos i/7r + x z 2ni Jg-i^ smr77r 

for some real 9 G (0, 1). From (|4.13p and (14. 15)) we can argue that 

(irfVa^)*'™ iLl v) /2L?\ Vt>0. (4.16) 

We notice that the equivalence in law (|4.16l) is independent of t as the formulae (|4.13p and (|4.15l) 

entail. The distribution h v o l v (x,t) of the process {r/L, t > 0} has Mellin transform (by making use 

L t 

of the formulae (|4.7[) and (|4.9|) ^ given by 

M[h u ol l/ (.,t)](r ) )=M[h v (;l)](r l )xM[l 1 ,(;t)} (^l + i)=i sm ^ t> f- 1 , t > 



f sm7r- 



with 5R{?/} € (0, 1). Thus, we can infer that 

~iu) law ~(v) i ~(v) 1 i~\ 



and h v o l„{x,t) — t l r(x/t) where r(w) is that in (|4.13|) . For the process i > 0} with 

distribution Z„ o h v (x, t) we obtain (from (|4.9[) and (|4.7|l ) 

X [|„ o fc„(., *)] fa) =X [!„(., 1)] (7?) X X [h v {; t)] (nU-V+l) = l ^W-W p-l 

v sin 707 



with ^{77} G (0, 1) and thus 
We have that l v o h v (x,t) = t~ 1 k(x/t) where k(x) is that in (|4.15|) 



L<# '^tx 1 L < p/ 2 I$\ t>0. 
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5 Main results 



In this section we consider compositions of processes whose governing equations are (generalized) frac- 
tional diffusion equations. When we consider compositions involving Markov processes and stable 
subordinators we still have Markov processes. Here we study Markov processes with random time 
which is the inverse of a stable subordinator. Such a process does not belong to the family of sta- 
ble subordinators (see (|4.12|) 'l and the resultant composition is not, in general, a Markov process. 
This somehow explains the effect of the fractional derivative appearing in the governing equation, see 
Mainardi et al. [19j. Hereafter, we exploit the Mellin convolution of generalized Gamma densities in 
order to write explicitly the solutions to fractional diffusion equations. We first present a new repre- 
sentation of the density law h v by means of the convolution e^ 1 introduced in Section [3] To do this 
we also introduce the time-stretching function </? m (s) — (s/m) m , m > 1, s € (0,oo). 

Lemma 2. The Mellin convolution e^Qs, (p n +i(t)) where fj,j — j v , for j — 1,2, ... ,n is the density 
law of a v-stable subordinator {r t , t > 0} with v = l/(n+ 1), neH. Thus, we have 

h v (x,t) = e* n (x,f n+1 (t)), x > 0, t > 0, v = l/(n+ 1), n e N. 



Proof From (|3.13|) we have that 



TT" r (1 — n + u-) 

M [ef{; Vn+1 {t))} fa) = fl^r^.) V n+iWV 1 - (5-1) 



From Gradshteyn and Ryzhik 12, formula 8.335.3] we deduce that 



k=l 

and formula (15.11) reduces to 



M [#•(., vn+i(*))] ( V ) = 3 (2n)n/2 ^ - (^(t))"' 1 . (5.3) 



Furthermore, by making use of the (product theorem) relation 

n-l / ,\ 

(see Gradshteyn and Ryzhik [l^, formula 3.335]) formula (|5.3[) becomes 

- [«■<-(«»] « = r( ^:;t;r- - 

(with 5i{?7} € (0, 1)) which coincides with (|4.7p . The claimed result is obtained. □ 

In light of the last result we are able to write explicitly the density law of a stable subordinator. 
For v = 1/2, Lemma [2] says that 

2.-1/2-1 e -*| 

h 1/2 (x,t) = e£(x,<p 2 (t)) = e 1/2 (x,(t/2) 2 ) = —^-^-^ , * > 0, * > (5.5) 
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which is the well-known density law of the 1/2-stable subordinator or the first-passage time of a 
standard Brownian motion trough the level t/y/2. For v = 1/3, from (|3.7p . we obtain 

1 i 3 / 2 / 2 i 3 / 2 \ 

ft 1/3 (a;,t) = (a;, <p a (t)) = e 1/3 *e 2/3 (x, (t/3) 3 ) = ^1 (^^TJ-y^J . £ >0,i > 0. 

For 1/ = 1/4, by (13.71) (and the commutativity under *, see Lemma [1}, we have 

h 1/4 (x, t) = ef{x, ipi(t)) = ex/4, * e 2/4 * e 3/4 (x, (t/4) A ) = e 1/2 * (ei /4 * e 3/4 )(x, (t/4) 4 ) 



where ^1/2(2) = a/ ir/2z exp{— z} (see [12t formula 8.469]). We notice that 

M [e* 3 (-,^ 4 (t))] (T))=M[h 1/2 oh 1/2 (;t)] (ry) 
which is in line with the well-known fact that 

#exp j-Aif ( ^ 2) j ^expj-A^f^} = cxp{-tA^ 2 }, 

< i>i < 1, i = 1,2. For v = 1/5, by exploiting twice (|3.7p (and the commutativity under ★), we can 
write 

h 1/5 (x, t) = ef(x, (t/5) 5 ) =(ei /5 * e 2/5 ) * (e 3/5 * e 4/5 )(x, (t/hf) (5.6) 

-STO^f ■- v - ljr »(VS jr »(^)* 

or equivalently 

^1/5(2:, *) = e* 4 (a;, (i/5) 5 ) =(e 1/s * e 3/5 ) * (e 2 / 5 * e 4/5 )(x, (i/5) 5 ) (5.7) 

= r/9 , (°° s- 1 '*- 1 !^ ( 2 x f^] K2 ( -L^ ds. 

55/2^2 3.2/5+1 y 5 V V^y J \5 5/2 

For ^ = l/(2n + 1), n e N, by using repeatedly (|3. 71) we arrive at 

u/2 t x/v-3/2 . . 

MM) = (*. M 1 ") , - > 0, t > 

where 



K,° v n {x,t) = I ...I K.„(x, si) . . .K v (s n -i,t) dsi . . .ds n -i 
Jo Jo 

is the integral (|2.5I) (as the symbol "on" denote) where n functions are involved and K. u (x,t) = 
x~ 2y ~ 1 K l/ (^y^t/x^j , x > 0, t > 0. We state a similar result for the density law l v and the convolution 

<?p'* ra (see Section[3]). Let us consider the time-stretching function ip m (s) — m s 1 /™, s <E (0, 00), m £ N, 
(ip = ip^ 1 where ip has been introduced in the previous Lemma). 

Lemma 3. The Mellin convolution g^ +1 ^* n (x,'ip n+ i(t)) where [ij — jv, j — 1,2,..., 71 and v = 
l/(n+ 1), neN, is the density law of a v-inverse process {l[ , t > 0}. Thus, we have 

l„(x,t) =4 n+1) '* n (z,V>n+i(*)), £>0, t>0, i/ = l/(n+l), neK 
Proof. The proof can be carried out as the proof of Lemma [2] □ 
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We obtain that li/^(x,t) = gf^ 2 (x, 2t 1 / 2 ) = e « /Vfti, x > 0, t > 0. Moreover, by making use of 
(|3~T1) and (15T21 . we have that 



V3^ = Ll*Kt (^ X — 



x > 0, t > 



il/3(^*)=iS/3*ffl/3(^3t 1/<, )- 7r Y V3 3/2 

and £1/4(2;, i) = * * 9i/a{ x ^ 4t 1 ^ 4 ) follows (thank to the commutativity under *) from 



33/4 * 52/4 * 9t/i(x, t) =g\ j2 * (p| /4 * 3i/ 4 )( a; > *) 



2 3 / 2 x 

7T t 



exp < — (sx) 4 



(sty 



ds 



where gi,^ 4 * gy 4 (x, t) is given by (j3Jj) and K 1 / 2 (z) = y/ir/2zexp{— z} (see 12], formula 8.469]). In 
more general setting, by making use of (|3.7I) we can write down 

fl V-.*(i/--i) (ljt) = _L_^^ Q°^(x,t), v = l/(2n+ 1), n 6 N (5.8 

where the symbol "on " stands for the integral (|2.5I) where n functions S i-» are involved and 



Q (x, *) = \2\J{x/t) l l v ^j , x > 0, f > 0. 



(5.9) 



Now, we present the main result of this paper concerning the explicit solutions to (generalized) 
fractional diffusion equations. We study a generalized problem which leads to fractional diffusion 
equations involving the adjoint operators of both Bessel and squared Bessel processes. Let us introduce 
the distribution ul'^ = g2 ° lv where g"^(x,t) = g^x^t 1 ^ 1 ) and the Mellin transform of which 
reads 

r f a=i + fi) r f a=i + 1) 

^RI'-ilK^— / i \ l-7M<3fiM< 1 + 7/^-7- (5-10) 

r( M )r(^ + i) 

We state the following result. 

Theorem 1. Let the previous setting prevail. For v = l/(2n +1), n G NU {0}, the solutions to 



D v t ul^ = Q 1 ^uY, x>0, t>0 
can be represented in terms of generalized Gamma convolution as 



(ttH 3 ") 1 ^ Jo 



Vv ( s , t)ds, x > 0, t > 



(5.11) 



(5.12) 



where Q-y u is the operator appearing in (|3.4[) . 



v u (s,t) 



I ... Qi^(s,si) . . .Qi^(s n ^i,t v /v)dsi . . .ds n ^i 
Jo Jo 2 2 



and Q i-i/ is i/iai in (|5.9p . Moreover, for v G (0, 1], we /icwe 



(1,1/); (Ai,0) 

(i,i); (M) 



(5.13) 



m terms of H Fox functions. 
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Proof. By exploiting the property (|2.6j) of the Mellin transform and the fact that 

x r >- 1 x f(x)dx = M[f(-)} ( v + 0), 



for the operator Q3.4p we have that 

M[g^ui>»(;t)} fa) 



(V - 7 + 1) + ~ - fa - 7) 

7 2 



=l(r? - - 7 )M [u^O, «)] (r? - 7) + - 1)(»? - 1)M [W,*)] fa - 7) 

=i(r? - 1)(77 - 1 +7M - l)M [u^O, t)] fa - 7 ) 
7 

where .M [^'"(-ji)] (77) is that in (15TTU)) . We obtain 



(5.14) 



M [Gy^i; t)] fa) = — fa - l)fa -l+W- 7)- 



(^+M)r(^ + i) ^ 



7 

-fa-1) , 
7 IW (a=i 



rfa)r 

r(^ + ,)r(^i) 



2t2=ii/ + l 

7 



r(H^ + M )r(fci + i 



r(M)r (a=i 



i/ - 1/+ 1 



[fl™ (■,*)] fa) 



and ul^ix^t) solves (|5.1ip for v £ (0, 1). In view of Lemma[3]we can write 

/•OC 

«2'' i (x J t)= / §2(x,a)^ / "'* (1/, '~ 1) (*,^(t))d8 



and by means of (|5 .8[) result (|5.12l) appears. Formula (|5 . 1 3[) follows directly from (|2.2p by considering 
formula (|3.10[) and the fact that 



jjra.TL 



1 



•' TTm,n 



(5.15) 

for all c G K (see Mathai and Saxena |2[j). □ 

We specialize the previous result by keeping in mind formula (|2.10[) and the operator (|3.4|) . 
Corollary 1. For v = 1, Theorem^ says that 



d 



1 f— a; 2 " 7 — - (7/1- 1)— a; 
\ 9a; 9a; 9a; 



^ 1- 



7 ,~.7,M 



tt 7 ' M , a; > 0, t > 0, 7 7^ 



dt 1 7 2 
where uj' 11 = g 7 is £/ie distribution of the GGP. 

This is because L\ a = t. Indeed, for v = 1, L\ is the elementary subordinator (see Q). 
Proof. If = 1, then the equation (|5.10[) takes the form 



T) — 1 



, R{t?} > 1 - 7/x 



(5.16) 
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where glix, t) = 17^(2;, i 1 / 7 ). For 7 > 0, we perform the time derivative of (|5.16[) and obtain 

I|-l/l)-7-l+7/l\ ^-7-1 \ 2=2=1 



7 V 7 / V 7 

=— (»? ~ tfiv - 7 - 1 + 7M)*t(»7 - 7) 

7 

which coincides with (|5 . 14[) and u7 ,Al (£, i) = <?7(:k, i), 7 > 0. Similar calculation must be done for 7 < 
and the proof is completed. □ 

Corollary 2. Let us write u^(x,t) = u];^(x,t). The distribution u^(x,t), x > 0, t > /1 > 0, 
f G (0, 1], solves the following fractional equation 

d v ( d 2 d \ 

w^ = { x w-^- 2 ^)<- (5 ' 17) 

In particular, for v — 1/2, we have 

x^' 1 f°° ( x s 2 1 
ufiJx, t) = —== / s M exp <^ }ds, x > 0, t > 0, u > 

which can be seen as the distribution of the process 

i G \B(2t)\> 1 > °} where B 

is a standard Brownian 

motion run at twice its usual speed and Gj'^ is a GGP. We notice that the process is a squared 
Bessel process starting from zero. 

Corollary 3. The distribution u 2 /^ — x > 0, t > 0, /i > 0, v G (0,1] solves the following 

fractional equation 

V 2 , u 1/9 2 9 (2 M -1) \ 2 , U 
In particular, for v = 1/3, we have 

M 1 ^(x,t)= / -TZTfiKi . ds, x>0, t>0, m>0 



7rr( M )v^7o «"- 1/2 3 V3 3/2 Vt 
and for /i = 1/2 we obtain 

^ 2{x ' t) = ^/tJ e ^ K i{w^l7i) ds> x>0 > t>0 

1/3 

which is the distribution of |-B(£ t )| where |-B(t)| is a folded Brownian motion with variance t/2. 

Acknowledgement The author is grateful to the anonymous referee for careful checks and com- 
ments. 
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